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The J/yz-hadron interaction is a key ingredient in analyzing the J/y/ suppression in hot hadronic 
matter as well as the propagation of J/i// in nuclei. As a first step to clarify the J/i//-hadron 
interactions at low energies, we have calculated i/\j/-n, J/y/-p and J/i//-nucleon scattering lengths 
by the quenched lattice QCD simulations with Wilson fermions for /3 = 6.2 on 24 3 x 48 and 
32 3 x 48 lattices. Using the Luscher's method to extract the scattering length from the simulations 
in a finite box, we find an attractive interaction in the S-wave channel for all three systems: Among 
others, the J/i//-nucleon interaction is most attractive. Possibility of the J/i//-nucleon bound state 
is also discussed. 
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1. Introduction 

A possible signature of the formation of the quark-gluon plasma (QGP) in relativistic heavy 
ion collisions is the suppression of J/y due to color Debye screening [jl|]. In the Pb-Pb collisions 
at CERN-SPS, the J/y suppression beyond the normal nuclear absorption has been discovered 
However, the data may be described either by the color Debye screening due to deconfined quarks 
and gluons or by absorption/dissociation due to comoving light hadrons. Recent BNL-RHIC data 
show a tendency that the J/y suppression is almost independent of collision energies between 62 
GeV and 200 GeV. The magnitude of the suppression is less than those predicted from color Debye 
screening or from the absorption by comovers, which may be understood by the recombination of 
c and c 

To understand the mechanism of J/\j/ suppression in those experiments, we need more precise 
knowledge on the J/i//-hadron interactions. In this report, we show our recent results of J/i//-hadron 
scattering lengths calculated in quenched lattice QCD simulations. (See [Q, ^ 0, [sj] and refer- 
ences therein for other approaches.) 



2. Formulation 

Suppose we have two hadrons in a finite box. The effect of their interaction appears as an 
energy shift AE relative to the non-interacting case. AE may be extracted from the correlator ratio 
R(t) for large t; 

R(t) = Gj /y; g(f) „ e -(E-my v -m H ); = -AE-t (2 1} 

where Gj/y. H (t) and G#(f) are the J/i//-hadron four-point function and the hadron two-point func- 
tion, respectively. 

In our calculation, we consider three scattering processes, J/y/"- 71, J/y-p and J/i//-./V(nucleon), 
which are most important for J/y absorption by comoving hadrons in relativistic heavy ion colli- 
sions. Since each hadron has spins, we need to make spin projection to good total spin states of 
J/i//-hadron two body system. For example, the J/y-N case reads 

G ]/v . N (t) = G xl2 {t)P 1 ' 2 + G 3/2 (t)pV 2 , (2.2) 

where G^nnm denotes the four-point function with a good spin quantum number (J = 1/2 or 3/2) 
and p'/ 2 ( 3 / 2 ) j s the spin projection operator 

The Liischer's formula tells us a relation between AE and the scattering observables such as 
the scattering length and the scattering phase shift. For the S-wave scattering phase shift, it reads 



[10] 



7T 3 / 2 fa 

tanfr(g) = - (2.3) 

with sign convention as which negative phase shift corresponds to repulsion. Here Zoo is a gener- 
alized zeta function defined by 
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p and L are the relative momentum of the two hadrons and the spatial size of the box, respectively. 
Here it is assumed that the interaction range is finite. Outside the interaction range, the total energy 
of the system is related to the relative momentum p as 



v/v 



m 2 + p 2 + a fm% + p 2 =E. 



(2.5) 



Here positive (negative) p 2 conesponds to the repulsion (attraction). If there are no interactions be- 
tween the hadrons, p takes discrete momentum in the finite box as p 2 = (2ti/L) 2 - n (n = 0, 1 , 2, • • •)• 
If there are interactions, p receives an extra contribution p mi and then a momentum squared divided 
by (2n/L) 2 is no longer an integer. 

The S-wave scattering length is denned as oq = lim p ^o tan (/?)//? and is related to the zeta 



function through Eq.([2.3|) as 



a 



(2.6) 



n=0 



2Z ()(1, 

where the scattering length is assigned to be negative (positive) for repulsion (weak attraction). The 



hadron scattering lengths based on the formulas Eqs.(2.3) and (2.6) has been extensively studied 
for %-% and N-N systems in Refs. [11, |12| ]. 



It is important here to discuss the asymptotic behavior of Eq.(2.6) for large L for the purpose 
of analyzing the system with attractive interactions. The large L expansion of the right hand side 
of Eq.Q at q ~ leads to @ 



AE 



2nao 
M res L 3 



l + ci[ r 



+ C 2 



5 ), 



(2.7) 



with c\ = —2.837297 and C2 = 6.375183. Let us try to solve Eq.(2.7) in terms of ao for given AE. 
In the case that AE > 0, both the expansion up to 0(L~ 4 ) and that up to 0(L~ 5 ) always have real 
and negative solutions. On the other hand, in the case that AE < 0, the expansion up to 0(L~ 4 ) 
gives no real solution for 



AE < 



K 



2\c] |M rPS L 2 ' 



(2.8) 



although the expansion up to 0{L~ 5 ) always has a real solution. 

This observation implies that some care must be taken to use the expansion especially for 
relatively strong attraction. Indeed, our lattice data show that the J/i// interaction with %, p and N 



are all attractive and the condition Eq.(|2.8[) is met for J/i//-p and J/\j/-N cases. Therefore, in our 
study, we use Eq.( |2.6| ) directly without the large L expansion to extract ciq. 



3. Results 



In our simulation, we employed unimproved Wilson gauge action and Wilson fermion. We 
have j8 = 6.2 on L 3 x T = 24 3 x 48 and 32 3 x 48 lattices with jc(charm) = 0.1360 and fc(light) 
= 0.1520, 0.1506, 0.1489. In the physical unit, the lattice sizes are L ~1.6, 2.1 fm, the lattice 
spacing is a ~0.067 fm, and ntjiy ~3.0 GeV and m % ~0.6-1.2 GeV. The number of quenched 
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Figure 1: The left (right) panel shows the quark mass dependence of the energy shift for the )ly-n (J/y-N) 
interaction in L = 32 lattice. The horizontal axis is the pion mass squared (m n a) 2 and the vertical axis is 
the energy shift AE. The open circles are the energy shift extracted from different K. The circles are the 
linear extrapolated points at the physical point. The open squares are the points which are evaluated with an 
assumption that energy shift is independent of quark mass. 



gauge configurations for smaller lattice is 161 and that for larger lattice is 169. Our error estimates 
are all based on the Jackknife method. 

In Figure 1, we show the quark mass dependence of the energy shift AE in the Jly-7i channel 
(the left panel) and the J/y-N channel (the right panel). The open circles are the energy shift 
extracted from the correlator ratio for different quark masses. The circles and the open squares are 
the results of a linear fit in quark mass and of a simple average over the data with different quark 
masses, respectively. The open squares are the estimate of AE without quark mass dependence as 
a reference. In both channels, we found that the interactions are attractive. 

In Figure 2, we show the volume dependence of the energy shift in the Jly-7i channel (the 
left panel) and the J/y-N channel (the right panel). The circles are the energy shifts from linear 
quark mass extrapolation to the physical point and the open squares are the results of using the 
constant quark mass dependence. Although the error bars are large in both channels, one can see 
the following tendency: In the Jly-7l channel, the absolute value of AE decreases as L increases. 
On the other hand, AE for J/y-N has opposite tendency. However, to make firm conclusions on 
this, we need to increase statistics and also collect the data for larger L. If it turns out to be true in 
high statistics data, one may conclude that J/y-n channel is attractive without a bound state, while 
J/y-N may have a bound state [p~3|]. 

Finally, in Figure 3, we show the volume dependence of the scattering lengths in the J/y-n 
(left panel) and the J/y-N (right panel) channels. The circles and the open squares are the results 
of linear quark mass extrapolation to the physical point and of fitting assumed constant quark mass 
dependence. To compare these absolute magnitudes of the J/i//-hadron scattering length with that 
in the 1 = 2 %-% channel, we put the empirical value of the %-% scattering length by crosses. Note 
that 1 = 2 71-n scattering is repulsive and we show its absolute value in the figure for comparison. 
The Jly-7i scattering length is negative and small compared to %-%. This is partly because the size 
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Figure 2: The left (right) panel shows the volume dependence of the energy shift from the J/y-7t(J/y- 
N) interaction. The horizontal axis is the spatial size L and the vertical axis is the energy shift AE. The 
circles indicate the energy shifts which are assumed to have linear quark mass dependence extrapolated to 
the physical point. The open squares show energy shifts estimated as if there is no quark mass dependence. 
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Figure 3: The left (right) panel shows the volume dependence of the scattering length in the J/1//-7T (J/i/z-AO 
channel. The horizontal axis is the spatial size L and the vertical axis is the S-wave scattering length ciq. The 
circles (open squares) indicate the scattering lengths which are assumed to have linear (constant) quark mass 
dependence. The crosses show the sign-flipped empirical value for the n-% scattering lengths in the 1 = 2 
channel as reference points. 
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of J/i// is small than the pion, and partly because only the gluonic exchange is allowed in the 
case. On the other hand, the scattering length for J/y-N could be order of magnitude larger than 
J/y-K, although the error bar is still quite large. 

4. Summary 

In summary, we study the J/i//-hadron scattering lengths by the quenched lattice QCD sim- 
ulations. We found attractive interactions in all J/y-n, J/y-p and J/y-N channels. Furthermore, 
the J/y-N scattering length is considerably larger than the J/i//-meson scattering length. Also, we 
found a sizable volume dependence of scattering lengths in all three channels. There is an opposite 
tendency of the volume dependence between J/y-n and J/y-N. There are several future problems 
to be examined further: To study whether the attractive J/y-N interaction could form a bound state, 
we need to have better statistics and simulations with larger lattice volumes, which is now under 
way. To confirm the validity of using the Liischer's formula, we need to check whether the potential 
range is small enough in comparison to the lattice size. Moreover, we need more careful analysis 
of inelastic contribution in our correlator ratio such as the D-D contribution in the J/\j/-p channel. 
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